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Abstract 

We consider the steady-state statistics of turbulence in the inertial interval. The Kolmogorov 
flux relation (4 / 5-law) is shown to be a particular case of the general relation on the current-density 
correlation function. Using that, we derive an analogous flux relation for compressible turbulence 
and a new exact relation for incompressible turbulence. 

PACS numbers: 
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I. INTRODUCTION 



The Kolmogorov flux relation is a rare exact analytic result one has for the correlation 
functions of the velocity v in the Navier-Stokes equation. It is traditionally derived consider- 
ing quadratic invariants like kinetic energy in an incompressible flow (or squared vorticity in 
two dimensions). Conservation of the kinetic energy, J v 2 /2, by an unforced Euler equation 
means that one can define the energy flux in fc-space and write the continuity equation for 
the energy spectral density [12] : 
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+ divll fc = 0, u(r) = v(r) - v(0). 



(1) 



Stationarity under the action of a large-scale force and small-scale viscosity means the con- 
stancy of the energy flux over intermediate wavenumbers [3, Q, 18]: 



ViiuiU 2 ) = -4Vi(^(r)^(r)^(0)) = -4e . 



(2) 



That can also be written for the third moment of the longitudinal velocity difference, ui = 
(u ■ r/r): 

(uf) = -12er/d(d+2) . (3) 

Here d is the space dimensionality. For d = 3 one obtains (uf) = — 4er/5, which is why the 
flux relation is sometimes called 4/5- law. Ana 
scalar turbulence, magnetized and helical flows 
the respective conservation law is non-quadratic (as the energy in compressible turbulence, 
for instance) then it is believed that the analogs of ([2]) are absent. In addition, there 
are no analogs of this relation for velocity moments of the orders n different from three. 



ogous relations are derived for the passive 
l3, 23, 25, 26, 27, 29|. If, however, 
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Experiments demonstrate that (u] 1 ) oc r^ n wit 
than one but smaller than n/3, see e.g. 1CJ, Il2 l. 



i the scaling exponents ( n , which are larger 



In this letter, we re-interpret the Kolmogorov relation in terms of currents and densities 
of the conserved quantities, which allows us to derive an analog for a compressible case and 
a new fundamental relation for an incompressible case. 
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II. GENERAL RELATION 



Consider a general class of classical field dynamics, 

d t q a + V-j a = f a , (4) 

where q a ,a = 1, N are densities, j a are currents and f a are the external random source 
fields. These equations describe local conservation laws and provide a canonical form for the 
effective hydrodynamic description of a slow evolution of large-scale perturbations. Since 
the zero wave-number component of the density is conserved, low wave-number components 



evolve slowly by continuity, see e.g. 



131 ] . The equations are closed via a constitutive relation 



that expresses currents in terms of the densities as a series in gradients: 

jt = m{p}) + E G tAip}) + ■ ■ ■ > ( 5 ) 

jb 

where dots stand for higher order terms involving more derivatives. The zero-order, reactive, 
term leads to a conservative dynamics while the first-order term describes dissipation. For 
fluid mechanics, consideration of higher order terms in Eq. §5§ is usually unnecessary so we 
limit ourselves to the following general class, which also contains the Navier-Stokes equation: 

w + = r - Wi (E ouip}) v jP b J . (6) 

We assume the standard mathematical formulation of the problem of turbulence where the 
orcing term f a is random and its statistics is stationary, spatially homogeneous and isotropic 



12j. This implies that the same properties hold for the steady state statistics of q a . The 
correlation length of the force will be denoted below by L. 

The derivation of the Kolmogorov type relation for Eq. (jH]) proceeds as follows. Consider 
the steady-state condition dt(q a (0,t)q a (r,t)) = (no summation over a). Employing the 
dynamical equation (jSj) and using statistical symmetries, one finds 

= d t {q a (0,t)q a (r,t)) = - 2 ^-{q a (0,t)F t a (r,t)) 

+2(q a (0,t)r(r,t)) - A (q a (0,t) ^ G^ h ({p}{r, t))V jP b (r, t) 

We consider the limit of large correlation length L of the forcing. That allows one to consider 
r much smaller than L yet still large enough, so that the last term becomes negligible as 



containing higher order of spatial derivatives. Because r <C L we have f a (r,t) ~ f a (0,t) 
and {q a (0,t)f a (r,t)) ~ (q a (0,t)f a (0,t)) = e a , where the constant e a is the mean input rate 
of {q a ) 2 - Hence we obtain for the single-time correlation function 

V,<g"(0)J?(r)) = e.. (7) 

Assuming in addition isotropy one finds 

(q a (0)i?(r)) = ^ • (8) 

III. PARTICULAR CASES 

A simple example is the passive scalar turbulence when some substance with the density 
q and diffusivity k is carried by a flow with the velocity v. The current is j = qv — nWq. 
The flux (Yaglom) relation is as follows: (q{0)q(r)v(r)) = er/d. 

Another example is the turbulence of a barotropic fluid where the pressure p(p) is a 
function of the fluid density p only. In this case, q = (pv, p), Fj = pViVj + p{p)Sij and 
j 4 = pv. The equations have the form 

dtp + V ■ (pv) = 0, d t (pvi) + dj(pViVj + pSij) (9) 
= -dj [Gi kb {{ P }) V kP b ] + f , (10) 

where we took into account that j 4 = pv is exact to all orders. The source often has the 
form f 1 = pVj$ that corresponds to an external potential $ (the analysis below can be 
easily generalized to other types of forcing as well). The exact form of kb is not important 
below but its presence is necessary for a steady state. Indeed, the rhs of Eq. ffTO]) breaks the 
energy conservation; the steady state holds due to the balance of the forcing that pumps 
fluctuations into the system and the dissipation that erases them. Now we use the general 
relation (jSJ) to derive the new relation for the compressible turbulence described by Eqs. (JS])- 
( TTO]) . The application of the relation to q A = p with e 4 = gives (p(0, t)p{r)v j(r)) = 
0. In fact this result holds for any relation between r and L. Indeed, the steady-state 
condition d t {p(0,t)p(r,t)) = = di(p(0)p(r)vi(r)) applied to the general (isotropic) form 
(p(0)p(r)vi(r)) = A(r)ri, gives A = Cr~ d , so that regularity at the origin requires C — 0. 
Application of Eq. (jHJ) to g J with j = 1,2,3 gives a non-trivial relation 

J2 (P(O)^(O) ^(rHrJ+pW^]) = ^ (11) 
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where e is defined in this (p(0)v(0) ■ /(O)) = e (we summed over j to get a more 

symmetric result). To the best of our knowledge, the formula ( II II) is new and presents 



a desired analog of Kolmogorov relation for the compressible turbulence see e. g. [21 ]. 
Probably, this relation was not derived before because it demands considering the steady- 
state condition for the fourth-order correlation function {p(Q,t)vj(Q,t)p(r,t)vj(r,i)), while 
usually in trying to find Kolmogorov type relations one considers steady state conditions for 



the second moment, like in magnetohydrodynamics by [25|. Note in passing that one can 
choose q as the energy density and obtain yet another relation analogous to (|T5|) below. 
In the incompressible limit, p = const and V ■ v = 0, the pressure term is zero; again, 



since it is a divergence-free vector that must be regular at the origin [see[22j. In this case, 
(fTTj) is reduced to ©, 

(v l (r)v J (r)v 3 (0)) = er i /d, (12) 

and Eq. ([3]) is implied. Hence (llip is indeed a general form of the Kolmogorov relation for 
an arbitrary Mach number. As we see, from a general viewpoint, the relation follows from 
the stationarity of the pair correlation function of the momentum density rather than from 
the energy spectral density. Indeed, e in (JTTj) is the input rate of the squared momentum 
and not that of the energy, which coincide (up to the factor 1/2) only in the incompressible 
case. 

We have shown how the Kolmogorov relation exploits the momentum conservation. Now 
in the same way we shall exploit the energy conservation and derive a new fundamental 
relation for incompressible turbulence. Energy conservation in the incompressible case means 
that the Navier-Stokes equation can be written as a continuity equation for the kinetic 
energy: 

d v 2 d 
di~2 ~ ~~c\i 

That means that one can choose q = v 2 which turns the general form (jSj) into a fifth-order 
relation. The only difference from the third-order relation (112j) is that the pressure term 
is now nonzero. Note first that the condition d(v 4 )/dt = gives the single-point pressure- 
velocity correlation function: 

(v\Vi P ) = (f ■ vv 2 ) + v {v 2 Vl S7 2 Vl ) . (14) 



v 2 



+ f-v + vvtfW (13) 
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For the different-point moment, Eq. (jSJ) gives 
d 



.. x v, ; \r) 



v 2 (r) 



+ p(r) 



^ 2 (0) ) = (f(r) ■ v(r)v 2 (0)) + v <t; 2 (0)^(r)V 2 ^(r)) . (15) 



This relation can be further simplified by decoupling small-scale and large-scale fields. Force 
and velocity are large-scale fields that change over the scale L so that one can put (f(r) ■ 
v(r)v 2 (0)) = (f-vv 2 ) at r C L as we did in deriving (jHJ). Velocity differences and derivatives 
are small-scale fields that change respectively over the scale r and the viscous length r\ see 
e. g. [3]. In particular, the local energy dissipation, z/t>j(r)V 2 fj(r), is a small-scale field 
with the correlation radius rj. When the distance r is in the inertial interval, r 3> 77, one can 
decouple large-scale and small-scale fields: v (v 2 (0)vj(r)V 2 f »(t*)) ~ —e(v 2 ). We also denote 
V = v(r) + v(0) and present 

4 (v t (r)v 2 (r)v 2 (0)) = - (u, [v 2 (r) - v 2 (0)} 2 ) = - ( Ui (u ■ V) 2 ) 

w - (u iU 2 ) (V 2 ) (cos 2 a) = -2{v 2 ) ( Ul u 2 ) = 8(v 2 )er/d . (16) 

Here a is the angle between u and V and in the second line we assumed that the vectors 
u and V are weakly correlated and can be decoupled when the distance r is in the inertial 
interval, r <C L. In the last equality we used the Kolmogorov relation ([2]). We thus find 
that in the inertial range, r\ <^ r <^ L, Eq. ( flBT) is reduced to 



A <^(r-Mr> 2 (0)) = (/ ■ vv 2 ) - 2e(v 2 ) = D . (17) 

We thus obtain the fundamental relation on the pressure-velocity correlation function which 
is the counterpart of the Kolmogorov relation ([3]): 

(v(r)p(r)v 2 (0)) = Dr/d. (18) 

That relation can be tested experimentally and numerically. From a formal viewpoint, 
(|11|17|) are particular cases of the Hopf equations which express the stationarity of the 
correlation functions. Hopf equations are generally not closed; they impose some relations 
between different structure functions but do not allow to find them see e. g. (3^]. On the 
contrary, our relations (11111171) allow one to find the correlation functions that are fluxes. 
Taking divergence of (fTBl gives (v 2 (0)t> j(r)Vip(r)) = D in the inertial interval, compare it 
with the single-point expression (|14|) . 
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J ^- (v 2 (r)v(0) [u 2 (r') - d [ti(r') ■ f'f J> = Dr/d. (20) 



Inverting the incompressibility condition, Ap = —'ViVj'VjVi, one expresses pressure: 

P[r) ~ (d-2)a d J \r-r'\ d - 2 ' ^ ~ F(d/2) ' ( } 

where Uj(r') = fj(r') — t>j(0). That allows one to present (fT8l) as an integral relation on the 
fifth-order three-point moment of the velocity field (we denote f = r'/r'): 

1 /" dr' 

We expressed it via the differences in the square brackets to make it explicit that the inte- 
grand is regular at the origin. 

Note that the force term in (fI7|) is the input rate of the squared energy which for a general 
force statistics cannot be expressed via the energy input rate. In the particular case of a white 
Gaussian force, such expression is possible: (fiViV 2 ) = (1 + 2/ d)e{v 2 ) and D = (2/d — l)e(v 2 ). 
Remark that the direct energy cascade, which is studied here, is absent for d = 2. For 
2d inverse energy cascade, one derives in a similar way: (v(r)p(r)v 2 (0)) = e^r/2 where 
€4 is the dissipation rate of the squared energy due to a large-scale sink. For the direct 
cascade of vorticity uo = V x v in 2d, analogous flux relations on the correlation functions 
(v(r)uj n (r)uj n (0)} do not contain pressure and were derived and analyzed by f|. 

Let us discuss the validity limits of (fl6|) where we neglected cumulants like the structure 
function (uu 4 ) c . Such structure function scales as r^ 5 and is subleading at sufficiently small r 
since £ 5 > 1. In other words, there exists the scale £, below which the decoupling is possible. 
We ask now how £ and (5 may depend on the only parameter that enters Navier-Stokes 
equation, the space dimensionality d. For d = 3 we expect £ ~ L. Usually, in statistical 
physics and passive scalar theory I^Q, Q the statistics is getting closer to Gaussian 
and decoupling improves when d grows. Somewhat counter-intuitively, one may expect a 
different behaviour in turbulence because of diminishing role of pressure. To understand 
how the pressure depends on d, consider how the identity Ap = — Vjf jVjVi behaves when d 
increases while the velocity components are kept fixed. The Laplacian has d terms of different 
signs. Assuming that in the limit d — ► 00 those terms can be considered independent, their 
sum grows like yd. The rhs contains d 2 terms and grows like d so that the pressure is 
expected to grow as p oc \fd. This is slower than v 2 oc d so that assuming d large one may 
neglect the pressure contribution into the correlation functions. Consider, for instance, the 
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tensor of the fourth moment of the velocity difference: 

(uiUjUkUi) = A(r)fifjf k ri + C(r) [SijS k i + 8 ik 5ji + 5 u 8 jk \ 
+B(r) [Sijfkfi + 5 ik r j r l + 5 u r k fj + Sjihri + bj k riri + 6^^) . 

At d — > oo the difference between longitudinal and transversal velocity components is ex- 
pected to decrease, as follows from the analysis of the second moment where incompressibility 
and isotropy requires (u\ — uf) = r(uf)'/(d — 1). Comparing (uf) = A + 6B + 3C with 
(u\) = 3C, we see that A, B <C C at d —>■ oo. Now, we generalize for d dimensions the 
relation between the velocity and pressure correlation functions, derived by [l^] for three 
dimensions, and remarkably find out that C-terms cancel out: 

That means that indeed (p 2 ) grows slower than (u 4 ) as d — > oo. Remind that generally 
the correlations functions containing p and u may scale differently yet considering the limit 
d — » oo at fixed r one may neglect the pressure term in (fT5l) . That would give 

V, <^(r> 2 (r> 2 (0)) = 2(f ■ vv 2 } - 2e(v 2 ) = Ae(v 2 )/d (22) 

(for Gaussian pumping) which is much different from e(v 2 ) given by ( Tl6l) . That means that 
at d — > oo there may exist an interval of scales, L 3> r ^> £, where cumulants are comparable 
to the reducible part so that i is a crossover scale between (fT7|) and ff22l . That would mean 
that i decreases with d while £5 — > 1 as d — > 00. Moreover, at d — > 00, the same analysis 
can be done for all odd moments: ( n = 1 for n > 1. We thus come to the suggestion that 
the scaling of incompressible turbulence in the limit d — ► 00 may be the same as the scaling 
of Burgers turbulence. This would not be that surprising since a single incompressibility 
condition imposed on d velocity components is expected to be less restrictive as d grows, so 
that flow configurations close to shocks give the main contribution to the moments. Note 
that for shock- like configurations, one cannot decouple u and V as done in deriving (IT7I) . 
Technically, neglecting pressure, makes all quantities J v 2n dr to be integrals of motion of the 
Euler equation so that the linear scaling of all odd velocity moments express the constancy 
of fluxes of these integrals of motion, like in Burgers, see e.g. Holder inequality then 
requires the linear scaling for even moments too, which corresponds to an extreme non- 
Gaussianity of the small-scale velocity statistics. Physically, pressure is a non-local field 
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which couples different regions in space and is expected to act like "an intermittency killer" 



as remarked by [20|, see also [5|, [24J. One way to find the exponents at finite d may be a 
large-c? expansion, which must thus be very different from that used by : one needs 

to start here from a Burgers-like limit rather than from a Gaussian statistics. That explains, 
in particular, why no substantial simplification was found and why pressure terms were not 
small in the large-c? perturbation theory that starts from a Gaussian statistics 11]. 

Needless to say that the tendency of the exponents to approach unity with d growing 
remains purely hypothetical on that level of analysis. However, if true, it means that the 
degree of non-Gaussianity of the statistics of a single velocity component grows with d, 
which agrees with the numerical comparison between three and four dimensions made in the 

n 

remarkable work by Gotoh and co-authors [6| . Scalar quantities made out of vector products 
contain the sum of d terms; one cannot generally conclude whether their statistics gets closer 
to Gaussian as d increases because of the competition between increasing non-Gaussianity 
of a single term and the averaging over d terms. 

To conclude, the formulas flllfl7lll8p is the main result of this work. 
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